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Abstract 

Classical  bifurcation  analysis  for  nonlinear  dynamics  is  used  to  derive  a  nonlinear  feedback  control 
law  that  eliminates  the  hysteresis  loop  associated  with  rotating  stall  and  extends  the  stable  operating 
range  in  axial  compressors.  The  proposed  control  system  employs  pressure  rise  as  output  measure¬ 
ment  and  throttle  position  as  actuating  signal  for  which  both  sensor  and  actuator  exist  in  the  current 
configuration  of  axial  compressors.  Thus,  our  results  provide  a  practical  solution  for  rotating  stall 
control  in  axial  compressors. 


1  Introduction 

An  axial  compressor  is  a  vital  part  of  turbine-based  aeroengines.  However,  the  engine  performance  is 
effectively  reduced  by  rotating  stall  and  surge  in  axial  compressors,  which  are  instabilities  that  arise  in 
the  unsteady  fluid  dynamics.  One  reason  that  these  unsteady  aerodynamic  instabilities  can  lead  to  large 
penalties  in  performance  is  that  they  are  difficult  to  predict  accurately  during  design.  Feedback  control 
has  to  be  employed  to  suppress  the  rotating  stall  and  surge  in  order  to  extend  the  stable  operating  range 
and/or  to  enlarge  domains  of  attraction  of  stable  equilibria  for  compressor  systems  and  to  improve  the 
engine  performance. 

As  such,  rotating  stall  and  surge  control  has  become  an  active  research  field,  especially  in  the  past 
decade  due  to  the  emergence  of  the  Moore- Greitzer  model  [12].  This  is  a  low  order  nonlinear  state- 
space  model  that  captures  the  nonlinear  dynamics  of  the  compressor  system  through  its  bifurcation 
characteristics  [2,  11].  Although  initial  research  has  focused  on  linear  control  for  rotating  stall  [5,  14], 
nonlinear  control  has  gradually  gained  more  recognition  and  attracted  more  attention.  In  fact,  the 
application  of  classical  nonHnear  dynamics  to  rotating  stall  and  surge  dynamics  motivated  a  simplified 
approach  to  rotating  stall  and  surge  control  based  on  bifurcation  theory.  This  idea  was  developed  by 

*This  work  was  supported  in  part  by  AFOSR  under  contract  no.  F49620-94-l-0415DEPSCoR,  and  by  WL/FIGC  of 
Wright- Patterson  AFB  under  a  Summer  Faculty  Fellowship  grant. 

^  Department  of  Electrical  and  Computer  Engineering,  Louisiana  State  University,  Baton  Rouge,  LA  70803-5901. 

^  Flight  Dynamics  Directorate,  Wright  Laboratory,  Wright-Patterson  Air  Force  Base,  OH,  45433-7531. 
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Abed  and  his  co-workers  [1,  2,  10,  16]  and  was  shown  to  be  effective  for  implementation  in  industrial 
turbomachinery  by  Nett  and  his  group  [6,  7]. 

In  this  paper  we  pursue  a  nonUnear  control  approach  based  on  bifurcation  theory.  The  proposed 
control  system  employs  pressure  rise  as  output  measurement  and  throttle  position  as  actuating  signal 
for  which  both  sensor  and  actuator  are  available  in  current  axial  compressors.  This  is  in  contrast  to  the 
linear  control  method  in  [5,  14]  that  requires  2D  (two-dimensional)  sensing  and  2D  actuation,  to  the 
nonlinear  control  method  in  [10,  16,  6,  7]  that  requires  measurement  of  amplitude  of  nonaxisymmetric 
disturbance  flow,  or  rotating  stall,  and  to  the  more  recent  work  in  [4,  8]  that  requires  air  jets  as  actuators. 
A  nonhnear  feedback  control  law  is  derived  based  on  classical  bifurcation  theory  for  nonlinear  dynamics 
and  is  shown  to  be  effective  in  eliminating  the  hysteresis  loop  associated  with  rotating  stall  and  in 
extending  the  stable  operating  range  of  the  compressor  characteristic.  In  comparison  with  the  existing 
results  on  nonlinear  control  of  rotating  stall,  our  feedback  control  law  accomplishes  the  same  goal  as 
those  reported  in  [10,  16],  and  more  importantly  has  the  potential  for  surge  control  as  well. 

The  schematic  diagram  of  the  axial  compressor  is  given  in  Figure  1.  The  notations  adopted  in  this 
paper  are  consistent  with  those  in  [12,  11]  with  minor  changes: 
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2  Moore- Greitzer  Model  and  Rotating  Stall  Control 


For  the  basic  compression  system  shown  in  Figure  1,  the  local  flow  coefficient  at  station  0  is  denoted  by 
6),  The  circumferential  mean  of  the  flow  coefficient  at  station  0  is  denoted  by  (l)(^)  and  is  given  by 

^(0=^^  =  '^(0  + 

where  with  the  flow  disturbance.  If  the  gas  flow  is  incompressible  and 

irrotational,  then  there  exists  a  disturbance  potential  flow  0,  t])  that  satisfies  Laplace’s  equation 


The  boundary  condition  is  taken  as  [11] 


dd"^  drf' 


=  0, 


where  If  ^h-  If  the  disturbance  flow  potential  and  circumferential  mean  of  the  flow  are  known,  then 
the  pressure  rise  coefficient  ^>(0  I®  determined  by  the  momentum  equation  [12]: 


1  {2^^  + 


9^  l„=o  2a  \  d^dr]  dOdrj 


The  pressure  rise  is  balanced  by  the  pressure  drop  in  plenum  governed  by  equation: 


“  </>r)  = 


rj, 


5  ^  =  Ft{4>t)  -  -Kt(I>t^ 


by  assuming  a  parabolic  throttle  characteristic.  For  uniform  and  steady  flow,  ^  the  performance 

characteristic  given  by  [12] 


=  H  C0  +  Ci(^-l)+C3(^- 


Laplace’s  equation  in  (1)  with  the  boundary  condition  in  (2)  yield  the  harmonic  form  for  disturbance 
potential  flow  [11,  1]  as  follows: 

CXD 

=  E  ,  ??  <  0.  (6) 

n=l 

Thus  the  flow  disturbance  at  station  0  has  the  form 

oo 

V’O  =  E  ^0”’  “  2VrA„(£)  sin(n^  -  r„£),  =  2n|a„|  smh{nlF)/W  (7) 

n=l 

with  Tn  a  constant  by  assuming  constant  speed  of  the  nth  harmonic  disturbance  flow,  i.e.,  constant 
wave  speed.  The  value  of  Ip  can  be  taken  as  infinity  provided  A„  has  a  finite  value.  Using  first  order 
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approximation,  or  single  harmonic  Galerkin  method  for  flow  disturbance,  Moore  and  Greitzer  obtained 
the  following  nonlinear  state-space  model  [12]: 
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for  compression  systems. 

If  the  following  changes  of  variables: 


$  = 


(8) 

(9) 

(10) 


(11) 


are  used,  and  with  notations  /?  =  2BH/W,  7  =  \/2H f  KtIW ,  and  a  =  1.5alc/{l  +  aXta.nh.(lp)/n),  then 
(8)  -  (10)  are  converted  into  the  following  simpler  form  of  ODEs  [11]: 

^  =  l($-$.^($)),  $.,($)  =  7^^  -  1,  7  >  0,  (12) 

$  =  -$-|-'P,($)-1-6C3$A2,  '$'e($)  =  Co  +  Ci$  +  C3$^,  (13) 

A  =  (tA  ^1  -  ,  (7  >  0,  (14) 


where  x  denotes  derivative  of  x  with  respect  to  “time”  variable  r.  The  throttle  position  7^  and  parameter 
(5  play  an  important  role  in  the  dynamics  of  axial  compressors.  The  cubic  polynomial  is  compressor 
dependent.  In,  [13],  the  following  values 


Co  =  8/3,  Cl  =  1.5,  C3  =  -0.5, 


(15) 


are  assumed,  while  in  reality  these  coefficients  are  uncertain.  It  should  be  clear  that  a  cubic  curve  is 
completely  parameterized  by  H ^  and  cq,  with  c\  =  — Scs  =  1/2,  that  define  its  shape  and  position. 
In  this  paper,  we  adopt  an  approach  that  the  variation  of  H  is  translated  into  those  for  cq,  ci,  and  C3. 
We  are  thus  led  to  the  following  assumption: 


1  <  Co  <  5,  Cl  >  0,  C3  <  0  and  Ci  +  3c3  =  0,  (16) 

for  the  coefiicients  of  the  performance  characteristic  curve.  The  parameter  a  is  less  important  and  less 
interesting. 

There  are  two  sets  of  equilibrium  points.  The  first  set  is  Ag  =  0  corresponding  to  the  axisymmetric 
flow.  In  this  case  the  compressor  performance  characteristic  is  given  by 

=  ^c(^e)  =  Co  +  Ci$e  +  Cs^l  (17) 

^Throttle  position  can  be  implemented  with  bleedvalves  if  used  as  actuator. 
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shown  by  solid  curve  on  the  upper-right  part  of  Figure  2  that  extends  to  the  left  with  dotted  line. 
The  steady  operating  point  of  the  compressor  is  determined  by  the  intersection  of  the  performance 
characteristic  'Fg  =  '®c(^e)  and  the  throttle  characteristic 

shown  with  two  dashed  lines  corresponding  to  two  different  7  values.  See  points  A  and  D  in  Figure 
2.  The  second  set  of  equilibrium  points  is  given  hy  Ae  ^  0  where  represents  the  amplitude  of 
the  nonaxisymmetric  disturbance  flow,  or  rotating  stall.  In  this  case,  the  performance  characteristic  is 
degenerated  into  the  rotating  stall  characteristic 


f  e  =  =  Co  +  {Ci  +  6C3)$e  “  5^3$^ 


(18) 


shown  by  solid  curve  on  the  lower-left  part  of  Figure  2  that  extends  with  a  sharp  slope  to  the  right 
indicated  by  dotted  line. 


Figure  2:  Schematic  of  compressor  characteristic 


By  employing  a  linearization  method,  it  can  be  easily  shown  that  the  solid  parts  of  ^ind 
are  stable,  but  dotted  parts  of  $5  are  unstable.  Thus,  if  one  gradually  decreases  the  throttle 

position  7  that  initially  determines  the  operating  point  D,  then  the  steady  operating  point  moves  along 
curve  to  point  A,  the  peak  of  the  compressor  performance  characteristic  curve.  Further  analysis 
with  center  manifold  theorem  for  nonlinear  dynamics  [15]  indicates  that  point  A  is  unstable  if  (15)  is 
used  for  in  (17).  Thus  nonaxisymmetric  disturbance  flow,  or  rotating  stall,  will  eventually  be 

born  and  the  steady  operating  condition  will  quickly  settle  at  point  B  on  the  rotating  stall  characteristic 
curve  $5.  There  is  a  substantial  drop  in  both  the  pressure  rise  and  flow  coefficient.  At  this  momemt, 
however,  increasing  the  throttle  position  does  not  recover  axisymmetric  flow.  Rather,  the  operating 
point  moves  along  curve  ^5,  resulting  in  even  lower  pressure  rise.  It  then  settles  quickly  at  point  D 
of  the  performance  characteristic  curve  after  reaching  point  C.  The  described  process  represented 
by  D-A-B-C-D  in  Figure  2  constitutes  the  hysteresis  loop  associated  with  rotating  stall.  Because  the 
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compressor  has  to  work  under  much  lower  pressure  rise  and  flow  coefficients  when  rotating  stall  occurs, 
prolonged  operation  under  this  condition  will  cause  severe  stress  and  overheating  in  rotor  blades,  and 
may  thus  damage  the  aeroengine.  However  little  can  be  done  to  return  to  axisymmetric  flow  once 
rotating  stall  is  fuUy  developed.  The  only  remedy  at  present  is  to  shut  down  the  engine  and  restart  it. 

Rotating  stall  in  axial  compressors  poses  a  great  challenge  to  both  turbomachinery  engineers  and 
control  engineers.  Because  of  the  existence  of  rotating  stall,  large  stall  margins  have  to  be  employed 
so  that  axial  compressors  can  only  operate  at  much  lower  pressure  rise  but  higher  flow  coefficient 
than  point  A,  resulting  in  poor  engine  performance;  because  of  the  existence  of  hysteresis  loop,  active 
control  may  not  be  effective  once  rotating  stall  is  fully  developed.  It  is  perhaps  beneficial  to  look  into 
the  dependence  of  equilibrium  point  (^ej^ej^e)  on  the  throttle  position  7,  and  to  understand  how 
local  stability  changes  with  respect  to  7.  Recall  that  the  throttle  position  synthesizes  the  effects  of 
combustion  chamber,  bleedvalves,  and  disturbances  to  the  flow  in  compressors.  Thus  it  is  a  highly 
uncertain  parameter  and  very  sensitive  to  inlet  disturbance. 

By  setting  the  right  hand  side  of  (12)  -  (14)  to  zero,  there  hold 

7^  =  =  0,  all  and  (19) 

Thus  equilibrium  solution  bifurcates  at 

'y  =  7c=  -7=^^==,  Ae  =  0,  =  Co  +  Cl  +  C3.  (21) 

That  is,  more  than  one  equilibrium  solution  is  born  at  7  =  7^.  Because  these  bifurcated  solutions  are 
independent  of  time,  they  are  called  stationary  bifurcation,  and  their  stability  can  be  determined  using 
the  well-known  linearization  method.  Clearly  rotating  stall  occurs  at  Ag  7^  0  that  is  also  a  function  of  7. 
Thus  stability  of  Ag  is  hinged  to  the  values  of  7  and  shown  by  the  bifurcation  diagram  (Ag,7)  as  in  (a) 
of  Figure  4  at  the  end  of  Section  4,  where  (15)  is  used  for  the  performance  characteristic  curve  as  in 
(17).  Again  by  convention,  the  solid  part  corresponds  to  stability  and  dotted  part  instability.  In  light 
of  [9],  this  is  exactly  the  subcritical  pitchfork  bifurcation  and  is  well  known  for  causing  the  hysteresis 
loop  in  classical  bifurcation  theory.  It  is  interesting  to  note  that  the  hysteresis  loop  persists  as  long  as 
Co  <  5,  but  disappears  for  cq  >  5  because  of  the  change  of  pitchfork  bifurcation  from  subcritical  into 
supercritical.  This  fact  was  pointed  out  by  McCaughan  [11].  A  more  interesting  observation  made  by 
Wang  and  Abed  [16]  for  the  case  of  cq  <  5  is  the  loss  of  stabilizability  with  linear  feedback  at  bifurcation 
point  7  =  7c.  Indeed,  setting  7  =  7^  +  and  linearizing  the  Moore-Greitzer  model  at  ($e,$e,Ae), 
Ag  =  0,  yield  the  linearized  model 
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where  —  $g,  X2  =  ^  sind  =  A.  Clearly  for  |$e|  ^  1?  the  unstable  eigenvalue  at  a(l  —  $g) 

is  not  affected  by  linear  control.  Thus  if  the  throttle  position  is  used  as  actuator,  linear  control  fails 
utterly.  This  observation  motivated  Liaw  and  Abed  [10],  and  Wang  and  Abed  [16]  to  use  the  nonlinear 
feedback  control  law  u  —  kA^  based  on  bifurcation  theory.  They  demonstrated  that  the  nonlinear 
feedback  gain  k  can  be  chosen  such  that  the  pitchfork  bifurcation  can  be  altered  from  subcritical  into 
supercritical,  thereby  stabilizing  the  bifurcation  point  j  —  and  >  0,  and  eliminating  the  hysteresis 
loop  associated  with  rotating  stall. 

While  the  nonlinear  control  law  in  [10,  16]  was  experimentally  validated  in  [3,  7],  it  remains  unclear 
how  A,  the  amplitude  of  nonaxisymmetric  disturbance  flow,  can  be  measured.  One  probably  has  to 
measure  local  flow  coefficient  (f)  around  the  circumference  of  the  compressor  in  order  to  estimate  A.  2D 
sensing  is  thus  required  especially  when  high  order  harmonics  are  involved.  At  present,  hot  wires  are 
used  for  measurement  of  local  flow  coefficient  [14].  These  are  very  delicate  devices,  making  it  difficult  for 
them  to  survive  the  hostile  environment  of  aeroengines.  Rather,  pressure  transducers  are  more  reliable 
and  endurable  in  the  volatile  flow  field.  In  addition  it  is  less  likely  for  the  control  law  u  =  kA^  to 
be  effective  for  surge  control,  especially  for  the  case  of  pure  surge  that  is  in  contrast  to  the  potential 
nonlinear  control  law  u  ~  A/'(^).  Therefore  a  challenging  problem  is  the  design  of  a  nonlinear  feedback 
control  system  where  pressure  rise  is  used  as  output  measurement  and  throttle  position  as  actuating 
signal.  It  should  be  clear  that  such  a  control  system  requires  only  ID  sensing  and  ID  actuation  even  if 
high  order  harmonics  are  involved.  Moreover  both  pressure  sensor  and  throttle  actuator  are  available  in 
current  configuration  of  axial  compressors.  The  difficulty  lies  in  the  design  of  nonlinear  feedback  control 
law  that  achieves  the  same  results  as  in  [10,  16,  3].  This  will  be  studied  in  next  section. 

3  Nonlinear  Feedback  Control  Based  on  Bifurcation  Theory 

This  section  focuses  on  the  design  of  a  nonlinear  feedback  control  system  where  pressure  transducer  is 
used  as  sensor  and  throttle  as  actuator.  As  discussed  earlier,  linear  control  fails  to  stabilize  the  equilibria 
of  Ae  >  0  for  the  case  cq  <  5.  Thus  nonlinear  feedback  control  law  has  to  be  employed.  We  begin 
our  discussion  with  the  geometric  shape  of  the  bifurcation  diagram  in  (a)  of  Figure  4.  The  bifurcated 
solution,  according  to  (19),  can  be  written  as  7  =  7(Ae)  for  Ag  >  0  at  least  locally.  The  hysteresis  loop 
associated  with  subcritical  pitchfork  bifurcation  is  due  to  the  two  maximum  values  of  7(Ae)  at  Ag  ^  0, 
whereas  supercritical  pitchfork  bifurcation  admits  only  one  maximum  value  for  7(Ag)  that  occurs  at 
Ag  =  0.  Intuitively  speaking,  the  elimination  of  the  hysteresis  loop  in  (a)  of  Figure  4  is  hinged  to  the 
elimination  of  the  two  extreme  values  of  7(Ae)  at  Ag  7^  0,  or  merging  them  into  a  single  one  located  at 
Ag  =  0.  The  nonlinear  feedback  control  law  proposed  in  [10,  16]  accomplishes  exactly  this  goal. 

For  the  control  system  considered  in  this  paper,  the  nonlinear  feedback  law  u  =  Af{^)  is  employed 
where  Af{*)  is  a  memory  less  nonlinear  function.  The  goal  is  to  design  such  that  it  eliminates 

the  hysteresis  loop  associated  with  rotating  stall  and  stabilizes  the  equilibria  of  Ag  7^  0.  Substituting 
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7  =  7o  +  A/’(W)  into  (12)  yields 


$  =  ^($+l-jV('5)v/¥-7o\/¥),  (23) 

where  7^  is  the  normal  value  of  the  throttle  position  that  may  not  be  the  same  as  7c,  the  critical 
value  at  bifurcation.  It  follows  that  the  nonbnear  feedback  control  law  u  =  M{^)  changes  the  throttle 
characteristic  curve  that  is  obtained  by  setting  the  right  hand  side  of  (23)  to  zero,  yielding 

(7c+A/'($e))"te  =  (l  +  ^e)^ 

In  general  the  steady  pressure  rise  \Pe  is  not  a  single- valued  function  of  $e  unless 

u  =  =  Kly/^  (24) 


is  chosen  with  K  the  feedback  gain  to  be  designed.  It  should  be  pointed  out  that  although  the  throttle 
characteristic  is  changed  for  K  ^  0,  the  compressor  performance  characteristic  $c  and  rotating  stall 
characteristic  remain  the  same.  We  will  derive  an  explicit  expression  of  K  that  eliminates  the 
hysteresis  loop  associated  with  rotating  stall,  and  stabilizes  the  equilibria  for  Ag  /  0  under  some  mild 
conditions.  For  this  purpose,  setting  the  right  hand  side  of  (13),  (14),  and  (23)  to  zero  gives  the 
expression 


($e  +  l-g)^ 

°  ^ci^e)  ’ 

(#e  +  l-g)^ 

°  ^.($e)  ’ 


for  Ae  =  0,  all  #e?  und 

(25) 

for  Ae  7^  0,  $e  =  ^1  -  A^. 

(26) 

These  two  equations  reduce  to  (19)  and  (20)  respectively  for  K  =  0  and  Jo  —  7-  We  will  first  demonstrate 
how  introduction  of  nonlinear  feedback  gain  K  as  in  (24)  can  actually  change  the  geometric  shape  of 
the  subcritical  pitchfork  bifurcation  in  (a)  of  Figure  4. 


Theorem  3.1  Suppose  the  coefficients  of^c  in  (17)  satisfy  (16),  and  cq  >  —4cs.  If  the  nonlinear 
feedback  gain  in  (24)  is  chosen  as 


_  2(co  IOC3) 
9C3  -  Cl 


then  Jo  —  7o(^e))  obtained  from  (26),  admits  only  one  extreme  value  for  —1  <  Ag  <  1,  and  it  takes 
place  at  Ae  =  0. 


Proof:  For  Ae  7^  0  as  in  (26),  direct  calculation  for  derivative  of  j^  with  respect  to  Ag  yields 


dj^ 


dAe 


27. 


djo  d$e 

d$e  dAe 


^7, 


Thus,  there  holds 


I'oi^e)  :  = 


djo  _  Ae  d7o  _  Ae 

dAe  ”  27o#e  d$e  ~  2$e'®'s($e) 


me), 


(28) 


8 


where  =  a/1  —  A^,  and  F($e)  is  given  by 


F($e)  =  2t/$e)  -  Ki^e)  ($e  +  1  “  K)  •  (29) 

It  is  claimed  that  F($e)  >  0  for  0  <  $e  <  1  provided  that  K  satisfies  (27).  Indeed,  there  holds 

=  ^si^e)  -  ($e  +  1  -  =  -dOcg^e  (a'  _  1  -  1$,  _  • 

Thus  for  C3  <  0  and  ci  +  6C3  =  3c3  <  0  by  Ci  =  — 3c3  as  assumed  in  (16), 


i$  I  (ci  +  6C3)  /cg  +  6c3 
2  30c3$e  “V  15c3 

for  aU  $e  >  0.  Moreover  the  gain  K  chosen  as  in  (27)  satisfies 


>  0 


j{  —  2(co  +  IOC3)  _  cq  +  IOC3  ^  ^ 


9c3  -  Cl 


6C3 


by  Co  >  -4c3  and  C3  <  0.  It  follows  that 


_  1  _  _ 


(ci  +  6C3) 

30C3$e 


<  A'  -  1  - 


I  Cl  +  6C3 


ir  =  K  -1 - 7=  <  0, 

15c3  x/5 


(30) 


using  the  hypothesis  ci  +  3c3  =  0.  By  again  C3  <  0,  the  above  inequality  leads  to  the  conclusion  that 
for  all  $e  >  0, 


(ci  +  6C3) 

30C3$e 


<  0. 


Evaluating  A($e)  at  =  1  gives,  using  (27), 


A($e  =  1)  =  2co  +  20c3  +  K{ci  -  9c3)  =  0. 

Since  /(#e)  =  A'($e)  <  0  for  all  >  0,  A($e)  is  a  strictly  decreasing  function  of  #6?  and  since 
A($e  =  1)  =  0,  we  conclude  that 

A($e)  =  F  >  0,  VAe  =  1/1  -  $2  >  0, 

by  $e  =  a/I  —  Al-  The  above  inequality  in  turn  implies  that 

I (  A  g c?'7o  f  ^  0,  for  <  0, 

^Ag  I  <  0,  for  Ag  >  0, 

by  7o  >  0,  >  0  according  to  (28).  Thus  7/Ae)  =  0  if  and  only  if  Ae  =  0  that  is  the  only  extreme 

point  for  the  function  7o(Ae),  in  fact,  a  maximum  point,  in  the  interval  |Ae|  <1.  I 

The  hypothesis  on  the  coefficients  of  the  performance  characteristic  curve  as  in  (16)  is  reasonable. 
Typically,  ci  =  1.5  and  C3  =  -0.5.  For  the  compressor  model  as  in  [13],  cq  =  8/3.  Thus,  the  feedback 
back  gain  in  (27)  is  obtained  as 

_  2(co  +  IOC3)  _  7 

9c3  -  Cl  9  ’ 
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as  required.  In  fact  by  examing  the  proof,  K  <  1  is  not  a  necessary  condition  for  the  result  in  Theorem 
3.1  to  hold.  Indeed,  (30)  indicates  that  K  <  1  can  be  improved  to  K  <  1  +  a/3/15  =  1.4472.  Because 
(30)  ensures  /($e)  <  0  for  aU  $e  >  0,  K  <  1.4472  is  again  not  a  necessary  condition.  This  wiU  be 
demonstrated  in  the  simulation  section.  However  because  large  K  implies  aggressive  control  strategy, 
and  consequently  causes  more  distortion  to  the  throttle  characteristic,  we  will  restrict  our  development 
to  the  scope  confined  by  0  <  iiT  <  1  in  the  rest  of  the  section. 

A  more  fundamental  question  on  the  nonhnear  feedback  gain  as  given  in  (27),  is  whether  or  not 
it  stabilizes  the  bifurcated  solution  at  jo  <  Ic,  and  Ag  >  0,  excluding  the  equihbria  Ag  =  0.  We  will 
consider  first  the  equilibrium  point  Ag  =  0  at  critical  value  of  7o. 


Theorem  3.2  Suppose  that  the  coefficients  of  the  performance  characteristic  curve  satisfy  (16),  and 
in  addition 

Co  >  -4c3,  C3  <  -1/2.  (31) 

If  the  nonlinear  feedback  gain  K  is  chosen  as 


K  =  max 


2  (cq  +  Cl  +  C3) , 


2(co+  IOC3)) 
9C3  -  Cl  / 


Then  0  <  K  <  1,  and  the  equilibrium  point  Ag  —  0  at 


is  asymptotically  stable. 


+  1  -  A  _  2-K 

\/t7  y/Co  +  Cl  +  C3 


(32) 


(33) 


Proof;  The  conditions  (16)  and  (31)  imply  that 


Co  +  Cl  +  C3  —  Co  —  2c3  >  —  6c3  >  3, 


2(co  +  IOC3)  —Co  —  IOC3  4c3  —  IOC3 


< 


9c3  -  Cl  -6c3  -6C3 

Thus  (16)  and  (31)  ensure  that  the  gain  chosen  as  in  (32)  satisfies  0  <  A  <  1.  Since  u  = 
nonlinear  feedback  system  has  the  form 


1.  (34) 

K/V^,  the 


Xi 

X2 

= 

.  ^3  _ 

'Vo 

-1 

0 


1(1- $2) 
0 


0 

0 

a(l-$2) 


Xi 

3^2 

.  ^3  . 

+  higher  order  terms. 


(35) 


where  ^2  =  ^  -  ^e?  and  =  A,  Because  the  equilibrium  solution  bifurcates  at  $e  = 

the  critical  value  of  7^  satisfies  (33).  Thus  center  manifold  theorem  [15]  can  be  applied  to  determine 
the  stability  of  Ag  =  0  at  70.  For  this  purpose,  set 


Xi 

Oil 

2  t 

'/3i  ' 

2^3  + 

J2  _ 

^2 

Oi2 

+  C>(a;|). 


Then  by  $e  =  1, 

X3  =  A  =  —a A?  —  2cr($  —  1)A  +  0(A*)  =  —ax^  -  2(tx2Xz  +  O(x^)  =  O(x^). 


(36) 
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Thus  there  holds 


2a2X^  +  3/322^3 


2aia:3  +  3/?ia:3  .  2aiX3  +  3/310^3 

_  ..o 


2020:3  +  3/322:3 


0{xl)  =  0{xf). 


On  the  other  hand,  the  nonhnear  feedback  system  in  (12)“(14),  in  conjunction  with  (24),  yield 


lo 


2/02^^  ^  ^ 
-1  0  X2 

/3“^  aixj  +  Pixl 

0  02x1  +  132x1 


+  0(0:3), 


by  (36).  Consequently  the  coefficients  of  lower  order  terms  0:3  for  A:  <  3  are  zero  that  results  in 


/3i  =  /?2  =  0,  «!  =  -3,  02  =  - 


2^/%' 


It  follows  that 


2:3  =  -o-a;3  -  2ax2X3  +  0{xf)  =  -<72:3  -  2ax3{a2xl)  +  O(x^)  =  -a{l  +  2a2)x\  +  0(x^).  (38) 

Direct  calculation  shows  that,  by  substituting  the  expression  for  K  as  in  (27), 


1  +  202 


-  6  +  3K  ^  2c3($e  -  1)  +  Co 
~  2C3 

2C3(^e  -  1)  +  Co  ^  -  1  -  Co  d  +  ^3  -  1  ^ 


due  to  C3  <  -1/2,  Cl  =  -3c3,  and  \Pe  =  ’^c(^e)  =  cq  +  Ci  +  C3.  For  K  =  2-  (cq  +  ci  +  C3)/3  =  2  - 
at  =  1? 

,  ,  ,  te  -  6  +  3/f  , 

1  +  202  = - - =  0. 

Hence  with  K  as  in  (32),  1  +  2o2  >  0.  There  are  two  different  cases  to  consider.  Case  1:  1  +  2o2  >  0. 
Then  the  local  asymptotic  stability  of  0:3  is  easily  concluded  by  the  center  manifold  theorem  in  light  of 
(38).  Case  2:  1  +  2o2  =  0.  In  this  case, 

-  ^  7o  _  1 

2’  ^/^■3 

by  (37).  Hence  in  order  to  determine  the  stability  of  ODE  in  (38),  higher  order  terms  are  needed.  Now 
(36)  gives 


xl  +  0{xl),  oi  =  -3,  02 


Similar  calculations  shows  that 


^  [v  ^  ~  6“^  ^  V  ^  ~  6^^ ) 


X2  =  -  (^ui  -  a:^  -  bixl  +  0(a;|)  =  0(a;|), 


X3  =  -a  [2a2  +  -jx3  +  0(3:3)  =  0(a:§). 
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Hence,  we  obtain 


«.  .  n  9  m  3  3  3 

6.  =  6,  =  0. 


It  follows  that 


„  1  5  3  5te  -  6  „ 

4  “  8  ~  4^e  ~  8^e  ^ 

by  $e  =  Co  +  Cl  +  C3  >  3  as  shown  in  (34).  Thus  the  ODE  in  (41)  is  asymptotically  stable,  that  in  turn 
implies  the  asymptotic  stability  of  Ag  =  0  by  center  manifold  theorem.  I 


The  condition  C3  <  —1/2  and  cq  +  ci  +  C3  >  3  as  in  (31)  is  true  for  the  compressor  used  in  [13].  In 
the  case  1  +  2cr2  =  0,  it  is  advised  to  choose  ^  >  0  such  that 

K  =  max  (2  -  ^  (co  +  Cl  +  C3) ,  ?(£o_+19£3)|  ^  ^  ^ 

I  3  '  "  9c3  -  Cl  J 

Such  a  choice  of  K  makes  1  +  2a2  >  0  that  validates  the  ODE  (38)  when  using  center  manifold  theorem, 
rather  than  (41).  This  is  due  to  the  fact  that  for  1  +  2a:2  =  0,  the  bifurcated  solution  jo  =  in 

the  neighborhood  of  Ag  =  0  has  its  derivatives 

^(0)  =  0,  i  =  1,2,3. 

Thus  the  curve  7^  =  7o(Ae)  is  very  flat  at  Ag  =  0.  Consequently  small  perturbation  of  7^  may  cause 
large  variation  in  Ag  that  is  not  desirable.  But  if  l  +  2a2  >  0,  then  7''(0)  7^  0,  minimizing  the  sensitivity 
of  Ag  with  respect  to  70  in  the  neighborhood  of  Ag  “  0. 

The  next  result  is  on  the  local  stability  of  Ag  7^  0,  but  close  to  Ag  =  0  when  feedback  gain  K  as  in 
(27)  is  used. 

Theorem  3.3  Suppose  that  the  coefficients  of  the  performance  characteristic  curve  satisfy  (16)^ 
and  (31).  If  the  nonlinear  feedback  gain  K  is  chosen  as  in  (32) ,  then  the  equilibrium  point  Ag  7^  0  is 
asymptotically  stable  provided  Ag  is  sufficiently  close  to  =  0. 

Proof:  Projection  method  as  in  [9]  can  be  used  to  show  the  local  stability  of  the  bifurcated  solution. 
We  shall,  however,  employ  a  more  elementary  proof  that  uses  no  more  than  linearization  method. 
Indeed,  linearization  with  xi  =  $  —  $e,  0:2  =  $  —  $g,  and  ojg  =  A  —  Ag  yields  approximate  linear  system: 


^2  “ 


3c3Ag  A  Cl  T  3c3  12c3Ag$g 


-2aAe$e 


-2(7  A^ 


The  corresponding  characteristic  equation  is  given  by 


+ 


2/32y^  +  (2(7  -  3C3)  Ag  -  (ci  +  3c3)J  a 


1  /  7o 

/32 


+  2(7 Ag )  (ci  +  3C3)  A 


12 


As  Ae  0,  all  the  coefficients  of  the  characteristic  equation  are  strictly  positive  except  the  constant  term 
that  goes  to  zero.  Thus  in  light  of  Routh-Hurwitz’s  stabihty  test,  the  local  stability  of  the  bifurcated 
solution  close  to  Ae  =  0  is  hinged  to  the  strict  positivity  of  the  constant  term.  Since  ci  +  3c3  =  0,  the 
constant  term  is 


/32 


30-70 

2/32^1; 


(8C3  -  lOcgA^)  Al 


> 


> 


/-j  67oC3\  2aAl  f  _  87^  \ 

2crAg(3A^  —  6)  +  (cq  +  Cl  +  C3)  ^  ^ 

Co  +  Cl  +  C3 


by  Cl  +  3c3  ~  0,  and  7^  <  (2  —  K)ly/Wl^  $e  <  cq  +  ci  +  C3  as  $e  <  1  for  Aq  ^  0.  This  concludes  the 
local  asymptotic  stability  of  the  bifurcated  solution  near  Ae  =  0.  I 


C3 


2  3  4  5 


Figure  3:  Allowable  region  for  coefficients  of  compressor  characteristic  curve 


Remark  1:  The  conditions  in  Theorems  3.1,  3.2,  and  3.3  amount  to  the  following: 

(i)  Co  <  5,  (ii)  C3  <  -0.5,  (iii)  cq  >  -4c3,  (42) 

by  Cl  +  3c3  =  0.  The  condition  cq  <  5  indicates  the  existence  of  hysteresis  loop  associated  with  rotating 
stall  for  axial  compressors  [11].  If  cq  >  5,  nonlinear  feedback  control  is  unnecessary.  The  conditions 
(ii)  and  (iii)  ensure  the  existence  of  the  nonhnear  feedback  law  as  in  (24),  such  that  0  <  AT  <  1  that 
changes  the  pitchfork  bifurcation  from  subcritical  into  supercritical,  and  stabilizes  the  rotating  stall 
characteristic  locally.  In  fact,  in  light  of  the  projection  method  as  in  ([9]),  our  nonlinear  feedback  law 
stabilizes  the  bifurcated  solution  Ae  7^  0  for  aU  0  <  $e  <  1-  This  is  due  to  the  fact  that  the  bifurcated 
solution  7o  =  7o(Ae)  does  not  admit  a  turning  point  for  $e  >  0.  The  conditions  in  (42)  are  represented 
by  shaded  area,  excluding  the  dashed  side,  in  Figure  3.  If  cq  and  C3  fall  into  the  shaded  area,  then  there 
exists  a  stabilizing  feedback  gain  K  that  can  be  computed  according  to  (32). 

Remark  2:  As  discussed  earlier,  A  <  1  can  be  improved  to  K  =  Amax  for  some  A  max  >  1  +  l/VS  for 
which  Theorem  3.1  still  holds.  If  local  stability  at  Ag  =  0  and  $e  =  1  is  the  only  concern  as  in  [10,  16], 
then  the  feedback  gain  in  (32)  can  be  replaced  by 

K  >  =  2  -  (43) 

by  the  proof  of  Theorem  3.2,  especially  (38).  In  this  case,  no  additional  constraint  is  required  on  the 
coefficients  of  compressor  performance  curve,  but  the  curve  7o(Ae)  defined  by  (26)  may  admit  more 
than  one  extreme  point.  Incidentally  for  the  compressor  in  [13],  Amin  =  7/9  as  well. 
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4  Numerical  Simulations 


The  nonlinear  feedback  control  law  as  presented  in  the  previous  section  is  tested  in  this  section  with 
numerical  simulations.  The  compressor  model  as  in  [13]  is  used  where  various  parameters  are: 

Az3l.75,  H  =  0A8,  B  =  2,  a  =  1/3.5, 

Co  =  8/3,  Cl  =  1.5,  C3  =  -0.5,  Ic  =  8,  Ip  —  oo. 

Figure  5  shows  four  simulation  plots  in  the  absence  of  the  feedback  control  law  (24).  Plot  (a)  is  a 
bifurcation  diagram  of  (Ae,7)  that  illustrates  the  subcritical  pitchfork  bifurcation  associated  with  the 
hysteresis  loop  in  rotating  stall  where  the  solid  line  indicates  stability  and  the  dotted  line  illustrates 
instability.  The  bifurcation  diagrams  in  (b)  -  (d)  of  Figure  4  are  obtained  from  (a)  of  Figure  4  using 
the  relations  satisfied  for  steady  equilibrium  solutions.  The  hysteresis  loop  in  (a)  of  Figure  4  clearly  has 
adverse  effects  for  using  the  throttle  as  control  actuator.  This  is  seen  from  the  fact  that  the  operating 
points  e?  ^e)  ^^e  not  single- valued  functions  of  7.  Thus  control  strategies  are  difficult  to  work  when 
rotating  stall  occurs.  However  after  using  the  nonlinear  feedback  control  law,  the  pitchfork  bifurcation 
in  (a)  of  Figure  4  is  changed  from  subcritical  into  supercritical  as  shown  in  (a)  of  Figure  5.  The 
bifurcation  diagrams  in  (b)  -  (d)  of  Figure  5  show  that  the  adverse  effects  resulting  from  hysteresis  loop 
are  eliminated,  and  ^e)  are  aU  single- valued  functions  of  7.  Hence  our  goal  stated  in  Section  1 

is  achieved. 

We  notice  from  (a)  and  (c)  of  Figure  5  that  at  the  critical  value  of  7,  small  perturbations  in  throttle 
position  may  cause  large  variations  in  the  amplitude  of  the  disturbance  fiow  and  the  pressure  rise.  This 
problem  can  be  eliminated  at  the  expense  of  using  large  gain  K.  Figure  6  shows  the  bifurcation  diagrams 
of  K  =  1.167  that  is  50%  increase  of  the  gain  K  =  7/9.  But  the  improvement  is  significant.  We  would 
like  to  emphasize  that  Theorem  3.1  holds  for  A"  >  1,  at  least  in  the  case  K  <  1  +  l/Vb  ^  1.4472 
as  discussed  in  the  previous  section.  Similarly,  the  stability  results  in  Theorems  3.2  and  3.3  hold  for 
K  >  Amin  where  Amin  is  given  in  (43).  The  effect  of  aggressive  control  gain  (A  =  1.167)  is  obvious  from 
(a)  and  (c)  of  Figure  6.  The  critical  value  of  7^  is  now  about  0.4,  35%  smaller  than  that  of  K  =  7/9, 
thus  increasing  the  tolerance  for  the  disturbance  in  throttle. 

5  Conclusion 

A  nonlinear  control  approach  based  on  bifurcation  theory  was  developed  to  design  a  feedback  control  law 
for  rotating  stall  control  in  axial  compressors.  The  proposed  control  system  employed  pressure  rise  as 
output  measurement  and  throttle  position  as  actuating  signal.  The  proposed  nonlinear  feedback  control 
law  was  shown  to  be  effective  in  eliminating  the  hysteresis  loop  associated  with  rotating  stall  and  in 
extending  the  stable  operating  range  of  the  compressor  characteristic.  In  comparison  with  the  existing 
results  on  nonlinear  control  of  rotating  stall,  our  feedback  control  law  accomplishes  the  same  goal  as 
those  reported  in  [10,  16].  A  more  interesting  problem  is  clearly  whether  or  not  the  same  feedback  law 
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works  for  surge  control  that  is  currently  under  investigation. 


Figure  4:  Bifurcation  diagrams  without  feedback  control 


Figure  5:  Bifurcation  diagrams  with  feedback  control 


Figure  6:  Bifurcation  diagrams  with  large  gain  feedback  control 
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